NMR relaxation rate and static spin susceptibility in graphene 
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The NMR relaxation rate and the static spin susceptibility in graphene are studied within a 
tight-binding description. At half filling, the NMR relaxation rate follows a power law as T 2 on 
the particle-hole symmetric side, while with a finite chemical potential fi and next-nearest neighbor 
t' , the (fi + 3t') 2 terms dominate at low excess charge 8. The static spin susceptibility is linearly 
dependent on temperature T at half filling when t' = 0, while with a finite fi and t', it should be 
dominated by + 3t') terms in low energy regime. These unusual phenomena are direct results 
of the low energy excitations of graphene, which behave as massless Dirac fermions. Furthermore, 
when S is high enough, there is a pronounced crossover which divides the temperature dependence of 
the NMR relaxation rate and the static spin susceptibility into two temperature regimes: the NMR 
relaxation rate and the static spin susceptibility increase dramatically as temperature increases in 
the low temperature regime, and after the crossover, both decrease as temperature increases at 
high temperatures. This crossover is due to the well-known logarithmic Van Hove singularity in the 
density of states, and its position dependence of temperature is sensitive to 5. 

PACS numbers: 81.05.Uw,71.10.-w,72.15.-v 
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I. INTRODUCTION 

Graphene, the latest carbon allotrope to be 
discovered^, is made out of carbon atoms organized 
into a honeycomb lattice. The characteristics of the 
honeycomb lattice make graphene a half filled system 
with a density of states (DOS) that vanishes linearly at 
the neutrality point, and an effective, low energy quasi- 
particle spectrum characterized by a dispersion which is 
linear in momentum close to the Fermi energy 2 -. These 
two features underlie the unconventional electronic 
properties of this material, whose low energy excitations 
behave as massless Dirac fermions^. 

There have been intensive theoretical and experimental 
studies on graphene to this date, for instance, half inte- 
ger and unconventional quantum Hall effect£&£'£, quan- 
tum minimum conductivity^^, bipolar supercurrentjii, 
ferromagnetismi 2 ., optical conductivity^ and the possi- 
bility of superconductivity^. On the other hand, the 
most interesting and promising properties from the tech- 
nological point of view arc its great crystalline qual- 
ity, high mobility and resilience to very high current 
densities^, the ability to tune the carrier density through 
a gate voltage^, as well as the fact that graphene exhibits 
both spin and valley degrees of freedom which might be 
harnessed in envisaged spintronic a 15 ' 16 ' 17 . For a review 
of other remarkable properties of such systems as well 
as a discussion of possible technological applications, the 
reader is referred to Ref.— Graphene is poised to be- 
come a new paradigm in solid state physics and materials 
science. 

Nuclear magnetic resonance (NMR) is usually an ex- 
cellent technique for probing the electronic properties of 
materials as it is sensitive to the DOS near the Fermi 
edge, and this method allows one to study static mag- 
netic correlations and low-energy spin excitations. For a 
material with a Fermi liquid state, the temperature de- 



pendent spin-lattice relaxation time T\ follows the well- 
known Korringa relation where 1/T\ varies linearly with 
temperature^. As one of the most powerful methods 
for investigating mechanisms of superconductivity of the 
many exotic types of superconductors being discovered 
today, it has turned out that the explanation for the 
peak in the NMR relaxation rate just in superconduct- 
ing state was that the DOS peaked dramatically at the 
edge of a energy gap 19 i 20 . Recent NMR experiments by 
Singer et al. showed a deviation from Fermi liquid behav- 
ior in carbon nanotubes with an energy gap evident at 
low temperatures 2 ^. In the framework of the Tomonaga- 
Luttinger liquid, the low temperature properties are gov- 
erned by a gapped relaxation due to a spin gap, which 
crosses over smoothly to the Luttingcr liquid behavior 
with increasing temperature^ 2 -. 

In spite of being few atoms thick, the system 
of graphene was found to be stable and ready for 
exploration^, and, it is believed that NMR should provide 
cornucopian and significant information on the electronic 
properties of graphene^. In the present paper, we study 
the NMR relaxation rate and the static spin susceptibil- 
ity of graphene within a tight-binding description. At 
half filling, the NMR relaxation rate follows a power law 
as T 2 on the particle-hole symmetric side, while away 
from half filling and with a finite next-nearest neighbor 
t', the (/i + 3i') 2 terms dominate at low excess charge 
6. The static spin susceptibility is linearly dependent on 
temperature T at half filling when t' = 0, while with a 
finite n and t', it should be dominated by (/x + 3t') terms 
in low energy regime. These unusual phenomena are di- 
rect results of the low energy excitations of graphene, 
which behave as massless Dirac fermions. Furthermore, 
when S is high enough, there is a pronounced crossover 
which divides the temperature dependence of the NMR 
relaxation rate and the static spin susceptibility into two 
temperature regimes: the NMR relaxation rate and the 
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static spin susceptibility increase dramatically as temper- 
ature increases in the low temperature regime, and after 
the crossover, both decrease as temperature increases at 
high temperatures. This crossover is due to the well- 
known logarithmic Van Hove singularity in the DOS, and 
its position dependence of temperature is sensitive to 5. 

The rest of the paper is organized as follows. The 
theoretical framework is introduced in section II. Our 
numerical result and discussion are shown in section III, 
and the paper is concluded with a summary in section 
IV. 



II. THEORETICAL FRAMEWORK 

Graphene is a two dimensional crystal of carbon atoms 
with a honeycomb lattice, which can be described in 
terms of two interpenetrating triangular sublattices, A 
and B, and then the electronic structure of graphene can 
be captured within a tight-binding descriptio n 18 i 23 ' 24 i 25 

H = -t^(al„bi +7la )+t'Y^(4<T a i+7v + b L b i+-Y<r)+ h - c - 
ir\<7 i-ya 

+ ^J2( a L a ™ + 6 LM' (!) 

ia 

where a,:.^ (aj CT ) annihilates (creates) electrons at the 
site Hi with spin a (a =t,l) on sublattice A, and bi t(T 
{b\ a ) annihilates (creates) electrons at the site R; with 
spin a (a =t, I) on sublattice B. t and t' are the nearest 
neighbor and next-nearest neighbor hopping energies re- 
spectively, and /i is the chemical potential. The presence 
of t' introduces an asymmetry between the valance and 
conduction bands, thus violating particle-hole symmetry. 
Specific values for t and t' have been estimated^ by com- 
paring a tight-binding description to first-principle calcu- 
lations. Following their estimates, we take t =2.7eV. To 
learn more on the effect of t' in graphene, cases with dif- 
ferent values of t' will be studied, and a typical t' =0.27 
eV will be paid more attention^ in the following. 

In the sublattice system, there are two coupled sublat- 
tices, and the energy spectrum has two branches. In this 
case, the one-particle Green's functions are matrices 

g(i-j,r) = ( 9 ^ i - j ' T \ 9a f.-i> T \), (2) 

where the longitudinal and transverse parts are defined 
as 

g mn (i - = -{T T mi(r)n](0)),m,ri = a,b; (3) 

with r is the imaginary time, and T T is the r order op- 
erator. Then the Green's functions are obtained as 

g aa {k,w) = - E ^-j^=9bb(k,uj), 

i/=l,2 w - £fc 

1 



respectively, where 
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= (4) 
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2[cos V3k x + 2 cos -k x cos (— -k y )]. 



From these, the DOS follows as 
p(u>) 
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with r — > + . For t' = an analytical expression for the 
DOS per unit cell can be derived^: 
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where Z(lj) — a /■f 1 with 
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and K(x) is the complete elliptic integral of the first kind. 
With a finite t', the DOS is evaluated by inserting unity 
in the form of an integral over the Dirac delta function 
as 



ku 

= ^£*(W(*) M ) 
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x 5 E / ^<5[e-(-l)V3+^]. (9) 
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Then, by interchanging the integration with k summa- 
tion, we obtain 

oo 

^) = ^E / de6[u-t'(gt-3)-(-iytg k }5[e-(-l) x tg k 

ku\_ OQ 

(10) 

with g k = \/3 + 7^. The second delta function enables 
us to replace g k by e/t(— 1) A , yielding to 

oo 
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where po(e) is the DOS with t' = 0, and then 
p(u>) = 



l + 4|;(w + 3t') „ * v^o 



c = |7[i+(-iryi+4^(^+3to 



(12) 



Now let us turn to evaluate the NMR relaxation rate 
and the static spin susceptibility. In general, 1/Ti mea- 
sures the local dynamics of the spins, and it is related to 
the transverse spin susceptibility xi(iw n ), which reads 
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X x(iWn) = - / dre^^(T T Sj(r)Sr(0)), (13) 



and after a straightforward calculation, we obtain 

(T T Sj(r)S-(Q)) = 9mn(0, r)g mn {Q, -r). (14) 
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FIG. 1: DOS for different values of the next nearest neigh- 
bor hopping i':(a) t'=0; (b) i'=0.27 eV; (c) £'=0.54 eV; (d) 
£'=0.81 eV with t=2.7 eV. 



With the help of spectral representation, we define 

A(k,Lu) = -2hn^2 9mn(k,uj), (15) 

m,n 

which has a relationship with the DOS 

^^A(k,u;)=^p(w). (16) 

fc 

In this representation, we can obtain the NMR relaxation 
rate finally 
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in which pg is the Bohr magneton, g is the electron g 
factor, h is the Planck constant, and j3 = l/hsT, where 
ks is the Boltzman constant. The F(k) are hypcrhnc 
form factors, which do not vary much with k in general. 
So we write all these form factors as F, which is indepen- 
dent of temperature. In the similar way, the static spin 
susceptibility can be derived as 

X = lim \\m x(q,u) 

q—>0 uj—>0 

g p B I p{u)~ i2 (18) 



4 cosh" (/3^) 



and x(q, w ) comes from 



X(i ~ 3,t) = J dre— <T r S?(r)S/(0)). (19) 
o 



III. RESULTS AND DISCUSSIONS 

NMR is a powerful method to characterize correlated 
states of materials as it is sensitive to the DOS near the 
Fermi edge. So we study the DOS firstly, and our results 
for different values of the next neighbor hoping t' have 
been shown in Fig.l. Focusing on the particle-hole sym- 
metric case, in Fig.l (a), it is clear that, besides the linear 
vanishing of the DOS at the Fermi level, there are marked 
van Hove singularities at the hopping energy, w = ±t. 
With a finite t', as shown in Fig.l (b), (c) and (d), these 
van Hove singularities shall appear at ui u = ±i, which is 
derived from Eq. 12, and we will show that t' plays an 
important role in graphene as it breaks the particle-hole 
symmetry. 

For both t' = and t' ^ cases, the DOS in graphene 
is markedly different from that in normal metals^ as 
its low energy excitations are two-dimensional massless 
Dirac fermions^, and the presence of t' shifts in energy 
the position of the Dirac point and breaks particle-hole 
symmetry. Our further results shall indicate that these 
kinds of important features of the DOS are at the origin 
of interesting properties of the NMR relaxation rate and 
the static spin susceptibility in graphene, as well as many 
transport anomalies in this material^. 

Having been familiar with the main features of the 
DOS in graphene, we now turn to the evaluation of the 
NMR relaxation rate and the static spin susceptibility. 
In Fig. 2 (a) and (b), the NMR relaxation rate and the 
static spin susceptibility as a function of t' for t=2.7 eV, 
T=10K at p=0 are plotted respectively. As shown by the 
dark line in Fig. 2 (a), the NMR relaxation rate follows a 
t' 2 power law at p = 0, while the static spin susceptibility 
is linearly dependent on t' as indicated by the dark line 
in Fig. 2(b). One of the most important properties of the 
NMR relaxation rate is its temperature dependent be- 
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FIG. 2: (a) The NMR relaxation rate and (b) the static spin 
susceptibility as function of t' for £=2.7 eV, T=10K at n=0. 
Dark lines indicate the exact numerical result and red lines 
indicate data computed from Eq. [21] or Eq. 1221 Inset: (a) 
The NMR relaxation rate and (b) the static spin susceptibility 
as function of temperature T for t=2.7 eV, t'=0, and /i=0. 
Dark lines with triangle indicate the exact numerical result 
and dash red lines indicate data computed from Eq. [21] or 
Eq. 1221 



havior. In the inset of Fig. 2(a) and (b), we have plotted 
the temperature dependence of the NMR relaxation rate 
and the static spin susceptibility respectively at t— 2.7 
eV, t'=0 and /i=0. As indicated by the dark lines with 
triangle, the NMR relaxation rate follows a T 2 power law, 
and the static spin susceptibility is linearly dependent on 
temperature T. 

With a finite chemical potential [i and t', the approxi- 
mate behavior of the NMR relaxation rate and the static 
spin susceptibility may be expressed analytically at low 
temperatures. From Eq |17i we see that only low energy 
part shall contribute to the NMR relaxation rate due 
to the properties of Fermi function at low temperatures. 
When \uj v \ <C t, as well as t' <C t, the DOS may be ex- 
pressed as 



2 a/3 \u} + 3t'\ 



3tt t 2 



approximately, and finally 
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FIG. 3: (Color online )(a) The NMR relaxation rate and (b) 
the static spin susceptibility as function of temperature for 
£=2.7 eV with £'=0.27 eV (red line), i'=0.54 eV (green line), 
£'=0.81 eV (blue line) at fi=0. Inset: Enlarge scale for £'=0.81 
eV. 



approximated by 



2 2 2y/3 1 
X-9 Ms^T - 7?i(^ 



3tt t 2 



3t') + 2fc s Tln[l 



-/3(/x+3t') 
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In this similar way, the static spin susceptibility may be 



The data computed within Eqf2T] and EqJ^U have also 
been shown in Fig. 2, which arc indicated by red lines, and 
each of them is very near to the exact result especially 
when t' < 0.2cV. 

The calculated temperature dependence of the NMR 
relaxation rate and the static spin susceptibility for dif- 
ferent values of t' at = are shown in Fig. 3 (a) and 
(b) respectively. At first glance, it seems that the NMR 
relaxation rate and the static spin susceptibility are al- 
most independent of temperature for a fixed t'. In the 
inset of Fig. 3(b), we enlarged the scale for the case of 
£'=0.81 eV, and the NMR relaxation rate increases as 
temperature increases, following a T 2 power law. How- 
ever, the enhancement with temperature is rather small 
comparing with the whole trend. 

This unusual phenomenon is a direct results of the low 
energy excitations of graphene, which behave as mass- 
less Dirac fcrmions. In low energy regime, the DOS in 
graphene is linear around the particle-hole symmetric fill- 
ing, and vanishes at the Dirac point, while the presence of 
£' shifts in energy the position of the Dirac point. Math- 
ematically it is clearly seen through Eqj20] Hence, with 
a finite chemical potential /i and £', we can describe our 
result within Eqf2T]and EqfJ^lat low temperatures rather 
well. In the particle-hole symmetric case, namely, £'=0, 
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FIG. 4: (Color online) (a) The NMR relaxation rate and (b) 
the static spin susceptibility as function of the excess charge 
for £=2.7 eV with t'=0 (dark line), £'=0.27 eV (red line), 
£'=0.54 eV (green line) and £'=0.81 eV (blue line) at T=10K. 
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FIG. 5: (Color online) (a) The NMR relaxation rate and (b) 
the static spin susceptibility as function of S for t=2.7 eV, 
t'=0.27 eV at different temperatures: T=10K (dark line), 
T=100K (red line), T=200K (green line) and T=300K (blue 
line). 



the NMR relaxation rate follows a T 2 power law. How- 
ever, the T 2 term in the NMR relaxation rate is negligible 
with respect to (// + 3£') for realistic values of /i and £'. 
On the static spin susceptibility, it is linearly dependent 
on T at half filling when £' = 0, while with a finite fi and 
t', it should be dominated by (/i + 3£') terms in low en- 
ergy regime, which may be described by Eg 1221 very well. 
Our results also show that £' plays an important role in 
graphenc since it breaks the particle-hole symmetry and 
is responsible for various effects observed experimentally. 

Arguably, one of the most interesting and promising 
properties from the technological point of view is the 
ability to tunc the carrier density in graphenc through 
a gate voltage^. Now let's turn to study the case with 
a finite excess charge 5, and the chemical potential /i is 
determined by 

f°° 1 1 

/ p(w)h=7 5 -}dcj = S. (23) 

J -co eP( u -ri + 1 2 J K J 

With the help of the gate voltage, one can control the 
density and type ( n or p ) of carriers varying their chem- 
ical potential^. The calculated S dependence of the NMR 
relaxation rate and the static spin susceptibility for dif- 
ferent values of £' have been shown in Fig. 4 (a) and (b) re- 
spectively. The NMR relaxation rate and the static spin 
susceptibility are linearly dependent on |<5| when \S\ < 0.1 
except £'=0.81 cV. While 6 is high enough, it is interest- 
ing to find that there arc several prominent peaks. In the 
particle-hole symmetric case, peaks appear around the 



center at <5=±0.157. With a finite £', these peaks shall 
appear around the center at (5=0.140, —0.144 for £'=0.27 
cV, (5=0.136, -0.147 for £'=0.54 cV, and 5=0.144, -0.260 
for £'=0.81 eV. 

These intriguing phenomenon may be predicted from 
the behavior of the DOS in graphene directly. When 
the chemical potential is located at the marked van Hove 
singularities, where the DOS peaks dramatically, hence 
peaks shall appear in the NMR relaxation rate and the 
static spin susceptibility^. To learn more on these peaks, 
the calculated <5 dependence of the NMR relaxation rate 
and the the static spin susceptibility with £=2.7 eV, 
£'=0.27 eV at different temperatures are plotted in Fig. 5 
(a) and (b) respectively. These peaks in the NMR relax- 
ation rate and the static spin susceptibility decrease as 
the temperature increases, however, these peaks are even 
pronounced at 300 K, and therefore, there peaks should 
produce a distinct effect on temperature dependence of 
the NMR relaxation rate and the static spin susceptibil- 
ity. 

In Fig. 6 (a) and (b), we plot the NMR relaxation rate 
and the static spin susceptibility as function of tempera- 
ture at different S for £=2.7 eV and £'=0.27 eV. Around 
(5=0.12, we can immediately separate the data into two 
temperature regimes by a crossover: the high tempera- 
ture regime and the low temperature regime. The NMR 
relaxation rate and the static spin susceptibility increase 
dramatically as temperature increases in the low tem- 
perature regime, and after the crossover, both of them 
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FIG. 6: (Color online) (a) The NMR relaxation rate and (b) 
the static spin susceptibility as function of temperature T for 
5=0.11 (dark line), 5=0.12 (red line), 5=0.121 (green line), 
and 5=0.122 (blue line), 5=0.14 (pink line), 5=0.158 (dash- 
dark-green line), and 5=0.16 (dash-dark-red line) at i=2.7 eV 
and t'=0.27 eV. 

decrease as temperature increases, and the position of 
temperature dependent crossover is sensitive to 5. 

For a fixed 8, the chemical potential varies slightly as 
temperatures varies, so when the chemical potential as- 
sociated with the temperature is located at the van Hove 
singularities, where the DOS peaks dramatically, there 
shall be peaks in both of the NMR relaxation rate and 
the static spin susceptibility, which may divide the tem- 
perature dependence of the NMR relaxation rate and the 
static spin susceptibility into two temperature regimes. 
The temperature at the crossover decreases as 8 increases 
when 8 < 0.140. As 8 tends toward 0.140, the NMR 
relaxation rate and the static spin susceptibility shall 
decrease as temperature increases almost in the whole 
temperature regime for current parameters, for instance, 
pink lines in Fig. 6 (a) and (b); and afterward, the tem- 
perature at the crossover increases as 5 increases when 
8 > 0.140 (dash lines in Fig. 6 (a) and (b) are respect 
to this illumination ). We only analyze the case when 
5 > at t' = 0.27 eV; because, for other values of if or 



when 8 < 0, the essential nature of peak's appearing is 
the same. 



IV. SUMMARY 

In summary, we have analyzed the NMR relaxation 
rate and the static spin susceptibility in graphene, which 
show a behavior that is not that of a normal metal. In low 
energy regime, the DOS in graphene is linear around the 
particle-hole symmetric filling, and vanishes at the Dirac 
point, while the presence of t' shifts in energy the posi- 
tion of the Dirac point. Hence, the NMR relaxation rate 
follows a power law as T 2 on the particle-hole symmetric 
side at half filling, while away from half filling and with 
a finite t' , the (/i + it') 2 terms dominate at low excess 
charge. The static spin susceptibility is linearly depen- 
dent on T at half filling when t' = 0, while with a finite 
(i and t' , it should be dominated by (fj. + 3t') terms in 
low energy regime. The next-nearest neighbor t' plays an 
important role in graphene as it breaks the particle-hole 
symmetry and is responsible for various effects observed 
experimentally. These unusual phenomena are direct re- 
sults of the low energy excitations of graphene. 

The NMR relaxation rate and the static spin suscep- 
tibility are linearly dependent on the excess charge \8\ 
when \8\ is small, while at high 8, there is a pronounced 
crossover which divides the temperature dependence of 
the NMR relaxation rate and the static spin susceptibil- 
ity into two temperature regimes: the NMR relaxation 
rate and the static spin susceptibility increase dramat- 
ically as temperature increases in the low temperature 
regime, and after the crossover, both decrease as temper- 
ature increases at high temperatures. This crossover is 
due to the well-known logarithmic Van Hove singularity 
in the DOS, and its position dependence of temperature 
is sensitive to 8. These properties show that graphene is 
a new class of materials with an unusual metallic state. 
Since the electronic density is easily controlled by a gate 
voltage, these phenomena can certainly be tested exper- 
imentally. 
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